Anderson Orthogonality Catastrophe in Disordered Systems 
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We study the Anderson orthogonality catastrophe (AOC) in finite conductors with diffusive 
disorder. The disorder averaged logarithm of Xi th e overlap between the ground states before 
and after adding a static impurity, is found to depend nonmonotonically on the disorder. In two 
dimensions (lnx -1 ) oc In 2 N in the weak disorder limit, thus showing a stronger dependence on the 
number of electrons N than in the canonical AOC. A very broad tail of the distribution of x> found 
numerically, is a signature of the importance of a few-level statistics at the Fermi energy. 



The Anderson orthogonality catastrophe is ubiquitous 
in physics, underlying such disparate phenomena as the 
zero-bias anomaly in disordered systems, the fractional 
Quantum Hall effect and the metal-insulator transition 
. The seminal paper of Anderson Q addressed a sys- 
tem of non-interacting electrons , confined to a finite vol- 
ume. No randomness (or complex geometry) were con- 
sidered. The introduction of a static impurity modifies 
the many-particle ground state (GS) of the original sys- 
tem, such that the overlap between the original and the 
modified GS, x — I (^1^') | 2 , is bounded (from above) by 
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Here \n) and \m') are the single-particle eigenstates of the 
original and the modified system respectively, with ener- 
gies £„,£,„< counted from the Fermi level; I is referred 
to as the Anderson integral. For a short range impurity 
whose potential has the strength A it was found H that 
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where N is the number of electrons: hence a "catastro- 
phe" as N approaches the thermodynamic limit. 

Since the AOC appears to be a rather general phe- 
nomenon, it is only natural to study it under more gen- 
eral circumstances. The importance of the inclusion of 
disorder is self-evident || |) . Not only may the disorder- 
averaged value of the overlap be different from its value in 
a pure system, but it can also induce mesoscopic fluctua- 
tions of X which are not necessarily narrowly distributed. 
Furthermore, one may consider longer range impurities 
as well: in 2DEG systems impurities are located off the 
conducting layer, implying softer impurity potentials. 

In this work we consider electronic systems in the pres- 
ence of quenched impurities that give rise to diffusive mo- 
tion of the electrons, with an additional impurity further 



introduced into the system. We address two main issues 
concerning the AOC caused by the additional impurity. 

(i) We consider the average value of the Anderson inte- 
gral as a function of N, of disorder, characterised by the 
electron transport mean free path, t, and of the range, b, 
of the added impurity parameterized by 
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where v is the single-particle density of states. We find 
that in the weak disorder regime and for pp 1 < b <C L in 
2d (L is the linear system size) and £ < b <C L in 3d Q| 
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Here go — 2ir 2 vD = 2tt 2 v vpt/ d is the conductance of a 
d-dimensional cube of size £ d . Thus, in 2d the decrease 
with N of the overlap Xi Eq. (]l|), is faster than a power 
law, showing a stronger N dependence than in the orig- 
inal AOC for clean systems. This is due to the diffusive 
character of electron motion which increases the proba- 
bility of return to the added impurity and thus the sen- 
sitivity of the modified GS to it. The overlap decreases 
faster when disorder increases, while in the strong disor- 
der regime (localization) it should practically disappear 
- hence a non-monotonic dependence of % on go- 
There exists also a ballistic contribution to the overlap, 
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For b < pp 1 it goes over to the "clean" result (||). For 
Pp 1 < b < I it dominates in 3d, while in 2d, for the max- 
imal N possible in the weak disorder regime (In N < go, 
see below), the diffusive contribution ([|) prevails. 

(ii) We observe the appearance of a wide distribution, 
p{x), of the GS overlaps. Numerically we find it to be 
anomalously broad even on a logarithmic scale. We ar- 
gue that the tails of the distribution are defined by the 
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level statistics of the last occupied and the first unoccu- 
pied levels. By projecting the problem onto a two-level 
Hilbert space and using random matrix theory, we derive 
this distribution non-perturbatively and show it to be in 
good agreement with the numerical results. The under- 
lying physics of this part is closely related to a recent 
study of the AOC in the context of parametrically varied 
random matrices ||. 

The logarithm of the overlap integral, X, can be exactly 
expressed in terms of Green's functions: 
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where g(e) = G R (e;r,r r ) — G A (e;r,r'). Scattering from 
all the other impurities is implicitly included in the for- 
mally exact retarded and advanced Green's functions, 
G R and G A . The presence of disorder does not have a 
significant impact on the mean value of the overlap X pro- 
vided that the additional impurity is point-like i.e. b < Ap 
in Eq. (^). Taking the limit 6 — + 0, i.e. u(r) — S(r), one 
reduces the overlap logarithm (||) to 
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Then, in the leading order in 1/go, the average (vv) 
is given by the disconnected contribution, (y) (y). Ap- 
proximating {v{e)) by v for |e| < E/2 (E is the band- 
width) and vanishing outside the band, one finds the av- 
erage of (^) equal to the Anderson expression, (X) = 
2n 2 X 2 \n(E / A) , where A = l/vL d is the mean single- 
particle level spacing. 

When the range of u(v) is finite, the mean value of 
X is more strongly affected by disorder. Now one needs 
to take into account both the disconnected (short-range), 
(.9*(e)> (ff(e'))) and connected (long-range), {g* (e)g(e')) c , 
contributions. The former is given by P]lCf| 
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where R = |R| = |r - r'|. Substituting Eq. (|) into Eq. 
(^), and performing the spatial integration and subse- 
quently the energy integration yields Eq. (||). 

The long-range contribution is obtained from Eq. (^) 
by considering the impurity average over the advanced- 
retarded parts of g*g yielding 
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Substituting the expression for u, Eq. (^|), performing 
the spatial integration, replacing the summation over 



q by the integration and finally performing the inte- 
gral over ui we obtain Eq. (^J), where we have used 
\n{L/b) 2 = IniV - ln[(fc F 6) 2 /47r] w In TV for b < L jn). 
The incipient AOC becomes negligibly small when the 
impurity is delocalized over the entire system. 

We have thus found that once the impurity range ex- 
ceeds pp 1 in d = 2 or £ ^> pp 1 in d — 3, the main contri- 
bution to the overlap integral is the long-range term (Q) . 
In d = 2, it gives rise to the suppression of \ which is 
stronger than a power law in N . Furthermore, the condi- 
tion that the magnitude of this effect exceeds the original 
AOC factor (0) is 
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At the same time we require weak localization (WL) cor- 
rections (neglected here) to be small, i.e. L < £, where 
^ is the localization length. In the presence of time re- 
versal (GOE) symmetry £ Q ~ £ exp go |Q , and these two 
inequalities are only marginally compatible. On the other 
hand, for a broken (e.g., due to an applied magnetic field) 
time reversal (GUE) symmetry ~ ^expg 2 , i.e. there 
exists a regime for which both Eq. ([l0]) and the require- 
ment on the smallness of WL corrections are satisfied. 

Arguably the most intriguing result concerning (I) is 
the dependence on disorder. It is clear from Eq. (|J) that 
in the weak disorder limit (and finite impurity range) the 
exponential suppression of the overlap is enhanced as g 
decreases. On the other hand one expects the AOC to 
vanish in the strongly localized limit. This implies a non- 
monotonic dependence on disorder. Such a behavior is, 
indeed, found numerically (see Fig. [TJ). 

The analysis outlined above concerns with the average 
value of X. As we show below, the distribution function 
p(x) is remarkably broadly shaped, so that the average 
alone does not provide us with the full information on 
the behavior of x- 

Our numerical study addressed a tight-binding Hamil- 
tonian in the presence of diagonal disorder: the site ener- 
gies assumed a box distribution of width W. The nearest- 
neighbor hopping, t, was chosen to be constant in mag- 
nitude, but with a phase factor when a magnetic field is 
applied. The underlying lattice was a two-dimensional 
square of area L 2 with periodic boundary conditions. 
For each realization of onsite energies the Hamiltonian 
is diagonalized exactly, both in the absence and in the 
presence of an additional impurity, parametrized accord- 
ing to Eq. (||). The calculation is performed for up to 
5 x 10 4 ,3.2 x 10 4 ,2 x 10 4 ,1.6 x 10 4 , 1.2 x 10 4 different 
realizations for the sample sizes L = 16, 20, 24, 28, 32, re- 
spectively, and several values of disorder. 

The broad distribution of X (normalized by its average 
value (X)) is presented in Fig. || for different values of L 
with filling N/L 2 = 1/4 and for W = 3t and W = 2t. 
Here we consider an impurity with A = 1 and 6 = 2. 
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The distribution we find is quite remarkable. Clearly it 
is even broader than log- normal |L3]]. The extended tail 
of the distribution corresponds to small values of x (large 
values of Xj (X)) and gives the appearance of power law 
behavior. The range of the tail shown in Fig. ^ (before 
the statistical scatter of data points washes it out) is at 
least 15 standard deviations. Furthermore the distribu- 
tion seems to be governed by a single scaling parameter 
X/(X), and does not depend on system size or on disor- 
der. The distributions exhibit extended tails and obey 
the single parameter scaling both for the GOE and the 
GUE symmetries. Upon close inspection it appears that 
the details of the tail differ between the GOE and the 
GUE scenarios, the occurrence of atypical events been 
suppressed in the latter (a "smaller tail"). This is ac- 
counted for below. 

The clear deviation of the distribution's tail from a log- 
normal shape suggests that large fluctuations in X are a 
manifestation of few-level physics. This is supported by 
a direct diagrammatic estimation of fluctuations. The 
leading contributions to (l 2 ) and to the higher cumu- 
lants of X come from the same diagrams as those for the 
local density of states (LdoS), as is clear from Eq. (||). 
The LDoS distribution by itself is log-normal [Q ; the ad- 
ditional energy denominator in Eq. (^|) leads to powers 
of (still) logarithmic divergence in all orders of perturba- 
tion theory for the cumulants of X higher than those for 
the LdoS cumulants. This indicates that the distribution 
of X is even wider than log-normal. Although it makes 
a direct diagrammatic (or renormalization group) analy- 
sis of the fluctuations of X hardly possible, this infrared 



divergence indirectly indicates the importance of just a 
few energy levels near the Fermi energy. 

The final confirmation of the few level picture is pro- 
vided by a scatter plot of x vs. the levels spacing S be- 
tween the iV+lst and A^th level for the quarter filled case, 
Fig. |. In these plots the differences (GOE vs. GUE) 
hinted to in the distribution plots become much clearer. 
While there is a finite (small) probability of finding small 
Anderson overlaps (x < 1/2) for the GOE case, there is 
almost no probability of finding such small overlaps in 
the GUE case. At the same time small values of x are 
strongly correlated with small level spacings, leading to 
suppression of the tail of p in the GUE case. 

Since the level spacing at the Fermi energy is strongly 
correlated with small values of the Anderson integral, it 
is reasonable to try to develop an analytic description 
for the tail of the distribution from the behavior of a 
two-level system representing the last occupied and first 
unoccupied levels in the system [|| . Prior to the introduc- 
tion of the impurity the levels have the energy —5/2 and 
5/2 and wave functions ipN(r) and tpN+i{r) respectively. 
After the introduction of an impurity XS(f = 0), the two 
level system is represented by the following Hamiltonian: 
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One can bound x from above by x', the contribution to 
the Anderson overlap due to the iV-th level (before and 
after the introduction of the impurity). This is readily 
done by diagonalizinig the above Hamiltonian: 



<5 + A(|vW0)| 2 H<M0)| 2 ) 
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We further assume that the distribution of the level 
spacing 6 is given by 

P(S/A) - / ^A) exp[-(V4)(<VA) 2 ], GOE 
[ / ' \(32/^ 2 (J/A) 2 )exp[-(4/^)( ( 5/A) 2 ], GUE 

and the wave function distributions are P(^ieip(0)) = 
P(3m </>(())) = (A/27r) 1 / 2 exp(-iV5ReV'(0) 2 /2) for the 
GUE case, while in the GOE case there are only real 
components of the wavefunction (N ~ (ppL) 2 ). The dis- 
tribution of x' is given then by 

p( X ) = J^dS J D^v +1 D^ S( X ' - X), (13) 

where which depends on all the integration variables, 
is given by Eq. (|l2"|), and Dtpw(N+i) denotes an integra- 
tion over the real and imaginary parts of ipN(N+i)(ty- 
The integration is performed numerically using a Monte 
Carlo method Jl5| , and the results for the tail region are 
shown in the inset of Fig. |2| As can be seen, deep in the 
tail region p{x) fits the distribution reasonably well. 



To summarize, we have found three main peculiarities 
in the physics of AOC in diffusive 2d systems. First, the 
average value of the overlap between the original and the 
modified ground states decreases with N faster than any 
power law. Second, it shows a non-monotonic depen- 
dence on disorder. Finally, the logarithm of the overlap 
shows remarkably wide fluctuations with disorder, with 
long tails which are mainly due to the contributions from 
only a few levels near the Fermi energy. 
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FIG. 2.: The distribution of X normalized by its average for 
different system sizes L at quarter filling N = L 2 /4. Open 
symbols correspond to W = 2t while filled ones to W — 3t. 
Left panel: No magnetic field (corresponding to GOE statis- 
tics). Right panel: Finite magnetic field (GUE statistics). 
Note that the larger values of 1/ (X) correspond to smaller 
values of x- Inset: The tail region for the GOE statistics at 
quarter filling, L = 16, W = 3t. The line corresponds to the 
two-level theory. 
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FIG. 1.: ((I) disorder/ In N) N (solid line) 

and {(T) disorder / In 2 N) n (dashed line) vs. the disorder W 
for an L = 32 system around quarter filling. 




FIG. 3.: Scatter plot of X vs. 5 for L = 16, W = 3t at 
quarter filling with GOE statistics. Inset: GUE statistics. 
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